Qubit Teleportation and Transfer across Antiferromagnetic Spin Chains 
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We explore the capability of spin- 1/2 chains to act as quantum channels for both teleportation 
and transfer of qubits. Exploiting the emergence of long-distance entanglement in low-dimensional 
systems [Phys. Rev. Lett. 96, 247206 (2006)], here we show how to obtain high communication 
fidelities between distant parties. An investigation of protocols of teleportation and state transfer 
is presented, in the realistic situation where temperature is included. Basing our setup on anti- 
ferromagnetic rotationally invariant systems, both protocols are represented by pure depolarizing 
channels. We propose a scheme where channel fidelity close to one can be achieved on very long 
chains at moderately small temperature. 

PACS numbers: 03.65.Ud, 03.67.Hk, 75.10.Pq 



Introduction. In order to accomplish the main tasks of 
Quantum Information, a sizable amount of entanglement 
is needed [1] . In addition, the particles that share entan- 
glement must be accessed individually for measurements 
and, quite importantly, they must be well separated in 
space. 

Recently it was shown [2J that in some spin models at 
zero temperature (i.e. in the ground state) a selected pair 
of distant sites A and B can be highly entangled. In some 
cases sites A and B may be taken infinitely far apart still 
retaining a high amount of entanglement, a situation that 
was termed long distance entanglement (LDE). An exam- 
ple of this situation is given by the end-sites of an open 
S =1/2 dimerized Heisenberg chain. Even for moderate 
values of the dimerization this effect is strong enough to 
develop non local correlations, i.e. entanglement, between 
the end-sites of an open chain of infinite length. 

The main aim of this Letter is to explore the actual 
feasibility of quantum teleportation and transfer across 
spin 1/2 chains that exhibit LDE. Having in mind realis- 
tic optical lattice implementations of spin chains 0], we 
consider the principal cause of decoherence which is given 
by the temperature. Using the same schemes proposed 
in Ref. we expect the entanglement between A and 
B to deteriorate when the temperature becomes of the 
order of the lowest excitation gap A. As this gap, which 
originates from the boundary conditions, typically van- 
ishes when the length L of the chain increases, we are led 
to explore the tradeoff between temperature and chain 
length. 

As will be clarified throughout this paper, antiferro- 
magnetic chains with global SU(2)-invariance have sev- 
eral advantages. Typically, in these systems rotational 
symmetry is never broken. As a consequence the two- 
particle reduced density matrix pab (obtained by tracing 
the total p over all the Hilbert space except sites A and 
B) maintains SU(2) invariance, i.e. it is a Werner state Q 



in the language of quantum information. Werner states 
are described by a single parameter which can be taken to 
be (a z A a z B ) p = Tr {p AB o- z A o- z B ) G [-1,1/3]. The interval 
( cf a <t b)p e [ — 1) — 1/3) corresponds to entangled pab- 

At T = the density matrix is p = \G)(G\, with \G) 
the ground state, while at finite temperature it is given 
by the canonical density operator p = Z~^e~° , with 
(3 = 1/T (in units of ks) and Z the normalization factor. 
At low temperatures we can approximate the thermal 
density matrix by retaining only the ground state and 
the first excited states. On quite general grounds [5( the 
ground state \G) is a total singlet, while the first excita- 
tions are given by a spin one triplet |m) labeled by the 
total magnetization: i§* ot = m = —1,0,1. Then at low 
temperatures we can write 
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where Eq is the ground state energy and A is the first 
excitation gap. Notice that this approximation correctly 
maintains rotational invariance. The thermal reduced 
density matrix pab (T) of A and B depends only on the 
following average value 
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((l\a A a%\l) + 2(l\a^%\l)) , (2) 



which has been written exploiting the SU(2) invariance 
The form ([2]) is particularly useful for numerical densit 
matrix renormalization-group (DMRG) simulations 
since it involves only the computation of the lowest-state 
correlation functions in the sectors m = and m = 1. 

In the situations analyzed in jij where LDE is present 
in the ground state, the S = 1 triplet state is localized 
near the sites A and B. As we will show below, the entan- 
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glement in pab (T) is maintained until T becomes com- 
parable with the gap A, when the triplet state becomes 
non-negligible. We are then led to prefer systems with a 
large gap A. Quite generally however, open systems with 
a finite bulk correlation length £ develop mid-gap levels 
scaling exponentially with the system size A ~ e~ L ^ 
Q. On the other hand, systems with a diverging corre- 
lation length give rise to an algebraic decay, A ~ L~ a . 
The generality of this conjecture - that establishes a rela- 
tion between bulk correlation length and the decay of the 
mid-gap - is a challenging question that deserves further 
studies. 

For the above-mentioned reasons we propose to use an 
open 5 = 1/2 Heisenberg chain with different interac- 
tions at the endpoints 



H, 



chain 



= H C + J f 



■ 02 + Sb • S 



L-l 



(3) 



as depicted in Fig. [T] (system ACB). In such a system 
there is strictly no LDE in the thermodynamic limit, but 
for finite size one can always choose J p /J small enough 
so as to have arbitrarily large entanglement between A 
and B in the ground state. Moreover we checked nu- 
merically that in this system the first gap A scales only 
algebraically with the size of the system L: A ~ L _a as 
can be seen in Fig. [2] Note the slow decay of the gaps 
due to the small value of a (see inset). 
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Figure 1: Model Hamiltonian considered for teleportation 
(joint measure between S and A) and for transfer (switch- 
ing on 7 at a given time). 

Teleportation. Entangled Werner (SU(2) invariant) 
states have several advantages when used as a resource 
for quantum informational devices. As far as teleporta- 
tion is concerned, one can show [|| that the standard tele- 
portation scheme @ is the best over all possible schemes 
at least in the region where a better-than-classical fi- 
delity is achieved. In the standard protocol an unknown 
state £ at site S (see Fig. [T]) is teleported to site B 
by making a joint Bell measurement on sites S and A 
and transmitting the result of the measurement j to B 
where a unitary transformation is applied. If A and B 
share a pure maximally entangled (SU(2) invariant) state 
\<P~)ab = (\U)ab - UT)ab)/V2, then the state £ is 
transferred to B exactly. In a realistic situation, exter- 
nal noise of any kind turns the pure state \?P~}ab into a 
non maximally entangled mixed state pab- 

Using this protocol with a Werner state as resource, 
the fidelity of teleportation does not depend on the the 
outcome j nor on the state to be teleported. By repeating 
the experiment many times with the same input state, 
the teleportation process is represented by a quantum 
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Figure 2: Finite size scaling behavior of the lowest gaps. In 
the inset is plotted the scaling exponent fitted with the law 
A = cL~ a . The data were obtained with a DMRG code using 
400-500 optimized states and three finite system sweeps. 



channel mapping input states £ at site S into teleported 
states A (£) at site B In this case, the teleportation 
channel is given precisely by a pure depolarizing channel: 



A(O=0£+(l-0) -1 



(4) 



The parameter i? which identifies the channel - some- 
times called shrinking factor - takes the simple form 
d = —{o~a°~b)- Obviously, A turns into an ideal chan- 
nel when § — 1, i.e. when pab is the singlet \i/>~)ab- 
The fidelity of teleportation is 



/ = Tr(£A(0) = 



1 + i9 



1 - 



A Bl 



that indeed does not depend on the state to teleport. 

For our class of systems, i? is given by Eq. ((2]). 
When the temperature is increased from zero, it even- 
tually reaches a value T* , above which the thermal 
state pab (T) becomes separable. This occurs when 
(<Ja°~b)t* = —1/3, that gives 



log 



2{l\o* A o%\l) 



1 



-{G\a A a* B \G} -1/3 



(5) 

Typical values in our scheme are obtained using the two 
qubit singlet and triplet pure states, for which we get 
T* = A/ log (3) w 0.9A. The gap A and the correlations 
appearing in Eq. |5]) can be calculated numerically as 
functions of L and J p . In Fig. [3] we plot the results, 
obtained from DMRG simulations for a chain of L = 
50 sites. In view of an optical lattice experiment, these 
curves could serve to locate the working point to achieve 
the maximal possible fidelity. 
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Figure 3: Fidelity of teleportation between end sites A and B 
as a function of temperature for the Heisenberg model. The 
curves refer to different values of the interaction J„. 



State transfer. As suggested by S. Bose [ll| open spin 
chains can be exploited for transferring quantum states 
from one end to the other end of the chain. Let a chain 
of length L be described by the Hamiltonian -f/chain- F° r 
times t < the chain is in its ground state or possibly 
in a state of thermal equilibrium Pchain = Z~ 1 e -,3ffchain . 
At a time t = a spin-spin interaction jSs ■ Sa between 
the sender S (that stores the pure state to be transferred 
|£}) and site A is switched on and let to evolve with the 
Hamiltonian, as depicted in Fig. [T] 



H = -ffchain +jS s ■ S A - 



(6) 



After a given optimal time t*, the initial state |£) gets 
transferred to site B with fidelity /. 

We stress here the importance of dealing with anti- 
ferromagnetic interactions. In this case, elementary ex- 
citations typically have relativistic linear dispersion for 
small momenta, i.e. lo (k) ~ v\k\ where v is the effec- 
tive speed of light. On the contrary in ferromagnetic 
systems, as the one originally proposed in [ll|, the dis- 
persion of elementary excitations is generally quadratic 
for small momenta. This fact leads to dispersive effects 
which limit the fidelity of transfer. 

From a quantum information perspective, one can eas- 
ily show that the state transfer protocol with SU(2) in- 
variant systems is precisely given by the depolarizing 
channel given by ([4]). The unique parameter specify- 
ing the channel is given in this case by i9 = (a z B (t)) p , 
p = |t) (Tl ® Pchain, where the time evolution is accord- 
ing to the total Hamiltonian The calculation of this 
quantity in a strongly correlated system is a non-trivial 



task. However, an approximation scheme is possible for 
the models where we observed LDE (or quasi-LDE). 

Although the spins on A and B do not interact directly, 
they experience an effective interaction mediated by the 
system C . Due to rotational invariance, the model © is 
effectively mapped, at every perturbative order, onto an 
SU(2)-symmetric Hamiltonian for the sites A and B 



He. 



JcsSa ■ Sb, 



(7) 



This approximation holds when the energy splitting J c g 
caused by H c g is smaller than the typical gaps in the un- 
perturbed Hamiltonian He- On the one hand, we know 
from conformal field theory [l^] that finite-size gaps in 
He scale as JlT 1 . On the other hand, J e g is noth- 
ing but the singlet-triplet gap A. We have numerically 
checked that A scales as JL~ a , in the system flchain, as 
can be seen from Fig. [21 The correct prefactor has the 
form <E>(J p /J). From perturbation theory we know that 
$(x) w x 2 for small x. This means that we can reliably 
approximate the model (J3j) with the effective Hamiltonian 
0, provided that $(J P /J) < L a ~ 1 , i.e. J p < JlP*-^/ 2 
when J p is small enough. 

Our scheme of approximation reduces the state trans- 
fer protocol to an effective three site problem where the 
time evolution is unitary by means of the Hamiltonian 
H = jSs ■ Sa + J c ff S a ■ Sb- The average is done with 
respect to the ensemble po = |£)(£|s ® Pab, where pab 
is the most general mixed state which preserve SU(2) 



invariance, i.e. 



Pab = 



9 

-a a ■ o\b, 



and g = (er^cr^) that includes decoherence effects from 
the environment C as well as the effect of temperature. 
The density matrix evolves as 



Pit) 



-itH 



{\0(t\s®PAB)e 



itH 



The fidelity of the transfer from site S to site B at a 
given time t is / (t) — Ttsab (p(t) After some 

calculations we get 



fit) 



= U 22 + 4 -9) (JSf + 7 2 ) ~ 7 Jeff (19 + 10 5 ) 

— 2 (1 + g) uj [u>- cos (tu>+ /2) + u + cos (ict>_/2)] 
+ 3 7 J cff (20-l)cos(wf)}, (8) 



where lu± = lo ± ( J cff + 7) and lo = yj J'* s - J e ff7 + 7 2 - 
The maximal possible interference (constructive and de- 
structive) is achieved when the frequencies are commen- 
surate each other i.e. for 7 = J e g. In this case the fidelity 
reduces to 

/ it) = ij[25 - 2g - 6 (1 + g) cos (J eff t/2) 
36 

+ (6(7 - 3) cos (J cS t) + 2(1+ g) cos (3 J eS t/2)]. 



The first maximum of this function is attained at a time 
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The value g = — 1 represents the ideal case where we 
have a pure singlet pab — \iP~)('>P~\ab at our disposal, 
with t* = tt/lu. In the non- ideal case, the time for best 
transfer gets only slightly shifted by a value which in the 
worst case (g = 1/3) is 1.448. The maximum fidelity is 
then 



%/3 (4 3 2 + Ag + 3) 3/2 + 24g 2 + 66g + 33 
48(l + <7) 2 



/*=/(**) 



= l--(g + l) + -(g + iy + [(g + if) . (10) 

As expected, the transfer is perfect for g = — 1. However, 
the transfer fidelity remains very high for all the possible 
values g e [—1,1/3]. The lowest possible value /* = 7/8 
is attained at g = (maximally mixed case). Anyway, we 
must restrict ourself to the situation where the approx- 
imation of unitary evolution is valid, i.e. g ~ — 1, and 
correspondingly the transmission fidelity is very close to 
1. As in the teleportation case, we have considered the 
transfer of a state using a Heisenberg chain playing the 
role of system C. In Fig. |4] the optimal transfer fidelity 
is plotted as a function of the chain length L at tempera- 
ture T = and T = 1CP 3 J, for some values of J p . In any 
case, we find a more-than-classical transmission fidelity 
even for chains of length 100 sites. For obtaining these 
results, the existence of entanglement between the two 
distant sites A and B is crucial . 

Now, let us draw an additional consideration. First, we 
note that t* cx = £"/(/$), with a < 1 (see Fig. HJ) 
On the other hand, our scheme is expected to be valid 
under the condition $(J p /J) < L" -1 which implies t* > 
L/J, consistently with the "flying" qubit picture where 
the information is carried by elementary spin excitations. 

Finally, we mention that the transfer protocol may be 
used also for sharing entanglement between distant par- 
ties [ill ] . In our situation, we already have entanglement 
between distant parties A and B, but we can ask how it 
may be further increased. 

The idea is to start having a maximally entangled sin- 
glet pi n — \ip~)(ip~\xs at sites S and at an extra neigh- 
boring site X completely decoupled from the rest. Then, 
we send the S part of the input state p m through the 
quantum channel described by our transmission proto- 
col. At a certain time t*, we obtain an outcome state 
living on the pair of sites X and B 
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Figure 4: Transfer fidelity at optimal time as a function of the 
chain length L. The curves refer to some values the coupling 
J p between the probes (A and B) and the Heisenberg chain. 
Results are reported at both zero and finite temperature. 



where p — 3(1 — #) /4 is the so called error probability. 
In this state, the concurrence between X and B is 



C( Po 



max [1 - 2p, 0] = max [3/* - 2, 0] , (11) 



while original state pab had a concurrence given by 
C (pab) = max [-3/20 - 1/2, 0]. Using /* from Eq. QE) 
it is possible to estimate that the concurrence is in- 



creased, i.e. C (p 



> 



C (pab), where the equality 



holds only when g = —1 (it is not possible to increase 
the entanglement of a singlet). The minimum value is 
achieved for the completely mixed case 5 = where the 
concurrence is C (p ut) = 5/8 = 0.625. 

Conclusions. We have given an explicit evidence that 
open antiferromagnetic Heisenberg chains may represent 
good quantum channels for teleportation and state trans- 
fer. This result relies mainly on the possibility to entan- 
gle the two end-spins (quasi-LDE) by choosing an ap- 
propriate coupling J p . We have shown that, despite the 
smallness of the lowest gap, high fidelities of both tele- 
portation and transfer may be achieved, with a tradeoff 
between temperature and chain length. It is tempting 
to speculate about the possibility of reproducing these 
effects in optical lattice environments [13j. 
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